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^ ■ Abstract 

x. 

D , The phenomenological characterstics of the electron and muon 

Q_ ( ' components of EAS with the size 10 5 < N e < 10 7 are obtained by 

the GAMMA installation of the ANI experiment at the Mt. Aragats 
in Armenia at the observation level 700 g.cm~ 2 . The experimental 
results are compared with other experiments and with the simulation 
^ ■ carried out using the CORSIKA code. 



1 ^Introduction 



The particle physics and astrophysics aspects of air shower studies are closely 
interlaced. The problems of sources and nature of the primary radiation and 
features of the high energy hadron interactions cannot be solved separately. 
The history of attempts to clarify the origin of the appreciable change in the 
index of the shower size spectrum (the knee), which was found more than 40 
years ago [|IJ and confirmed by other experiments, is the best manifestation of 
that. Many new experiments were born devoted to the observation of various 
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components of EAS, including Cherenkov light in the atmosphere and experi- 
ments under water. However, up to now there is no unequivocal explanation 
of the reason of this phenomena. The most popular interpretation of the 
knee is the existence of the knee in the primary cosmic ray energy spectrum 
at energy about (3-5) 10 6 GeV and, as a possible consequence, the change 
of the mass composition. There are different experimental analysis leading 
to opposite conclusions about primary composition after knee (lighter ||, 
normal §§ or heavier §110). 

At the same time there is an alternative explanation of the knee, also 
based on some inconsistencies in the EAS data, connected with a modifica- 
tion in the hadronic interaction properties || to the change of the number 
and energy spectrum of the secondaries. The complex installations, which 
measure electromagnetic component as well as muon and hadron compo- 
nents of EAS have the best opportunities to explore the nature of the knee. 
HAD RON- 2 |g on Tien-Shan and KASCADE |If are installations of such 
a kind. The hadron component data obtained from these installations show 
the fast absorption of the high energy hadrons, which is impossible to be 
explained by the energy increase of the interaction cross-section only. At the 
same time, the electromagnetic component data are in an agreement with 
the QGSJET model and correspond to the mixed composition of primary 
with progressive changing of the mass composition. In the work |5| on the 
bases of the Tien-Shan data analysis it is supposed that there is no knee in 
the primary spectrum and knee in the EAS size spectrum is explained by 
the change of the hadronic interaction in the atmosphere. This explanation 
is based on the results of the EAS electromagnetic component study |12| , 



according to which there is no knee neither in the young showers spectrum 
(showers with S < 0.75, generated by primary protons interacted in the depth 
of atmosphere), nor in the old showers (S > 1.05) one. And the existence of 
the knee in the same spectra of other experiments could be explained by the 
missing of the young showers. 

On the other hand it is not so clear situation with the so-called reverse 
knee in the EAS size spectrum at N e ~ (2 — 3) * 10 7 , which was discovered 
in [13|]. Subsequently the same irregularity was obtained once again in Tien- 
Shan experiment [T2"| ]. 

At present, practically all experiments on research of the cosmic rays in 
the energy range 10 14 -10 17 eV are pointed to investigate the nature of the 



knee. We have to point out that the KASCADE JTTJ, CASA-MIA (TJ and 



EAS-TOP 1 15] experiments present the basic advantage to have their data 
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analyzed using multiparameter procedures. Only such an approach to study 
the EAS characteristics allows some advance in understanding of the reasons 
of the existence of the knee in the spectrum on N e . In this paper, we present 
the experimental results of the GAMMA installation of the experimental 
complex ANI also working in this field. 

The GAMMA installation is located on hillside of the Mt. Aragats in Ar- 
menia (3200 m a.s.L). In comparison with the other large scale experiments 
with the same goals, the high-altitude observation level of the GAMMA in- 
stallation provides some advantages, in particular for minimizing the intrinsic 
fluctuations of the observables due to the stochastic character of the EAS de- 
velopment in atmosphere. This will add some new observations, which would 
clearly complement the data from KASCADE (near sea- level), CASA-MIA 
(1200 meters), and EAS-TOP (2000 meters). In the present work, show- 
ers are simulated using the CORSIKA code 5.20 [TBI in which the hadronic 



interactions are described by the VENUS model, [0. On the other hand, 
to compare the simulation data with experimental ones, the normal mixed 
composition was used, i.e. proton-40 %, a-21%, light-nuclei (<A>=14)-14%, 
medium-nuclei (<A>=26)-13% and heavy-nuclei (<A>=56)-12% 



Present status of the GAMMA EAS instal- 
lation 



m 



The GAMMA installation |18| was realized as a part of the project ANI [[19 
an attempt to continue experimental study of the hadronic interactions and 
the energy spectrum and mass composition of the primary cosmic radiation 
in the energy range of 10 14 — 10 17 eV. 

After some years spent to enlarge the effective area of the muon un- 
derground detector, which was necessary to elaborate methodical studies of 
the detector parameters, to investigate carefully the array response and to 
determine the precision of the shower parameter estimation, the GAMMA 
experiment is, now, effectively running. 

GAMMA is a central type array and consists of two main parts (figure 1): 

(i) the surface part, for the registration of the EAS electromagnetic com- 
ponent; 

(ii) the large muon underground detector, to register the EAS muon com- 
ponent. 
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2.0.1 The surface scintillation array 



The surface scintillation array constists of 25 groups of 3 scintillation detec- 
tors placed on concentric circles with radii of 17, 28, 50 and 70 m. Each of 
the detectors has 1 m 2 of effective area. They are distributed on the full area 
of «1.5-10 4 m 2 . An additional station with 20 detectors of the same type is 
placed at the 135 m from the installation centre. The energy threshold for 
the registration of the electrons flux, E l ^ res ~ 9.5 MeV, is estimated taking 
into account the thickness and the specificities of the detectors and registra- 
tion stations. The calibration of detectors has been carried out by periodical 
(each hour) registration of the cosmic ray background particles . It was taken 
into account that the most probable value of the registered particles differs 
from 0° and depends on the observation level |2(| . 

Each of the 25 registration stations is equipped with the timing channel 
which allows to determine the angular coordinates of the shower axis. The 
trigger condition for one EAS registration is that in each of four groups of 
detectors placed at 17 meters from the array centre, two detectors have to 
register a flux density with p > 3 part/m 2 . 

The main EAS parameters: coordinates and angles of the shower axis 
(A, Y, 9, <p), number of electrons (N e ) and age parameter Snkg are obtained 
using the SPACE code [31]] created on the basis of the statistical methods 



for the solution of the inverse problem |22fl . Taking into account the detector 
response, GAMMA array geometry and using Monte-Carlo procedure the 
data bank of pseudo experimental showers was created and treated by the 
same SPACE code. 

As the result of this study, the following accuracies on the EAS parameter 
estimations are obtained: AA, Ay ^ 3m (for R < 40m), o$ ±2 1.5°, ^ 8°, 
a N J < N e >< 20%. 

Furthermore, the analysis of the shower registration efficiency shows that 
the showers with particle number N e > 3 ■ 10 5 within 20 m and N e > 10 6 
within 40 m from the array centre are registered with 100% efficiency. 



2.0.2 Underground muon detector 

The layout of the muon underground detector is shown in figure2. As for 
the GAMMA surface part identical scintillation detectors with 5* = lm 2 
were used for the registration of the EAS muon component. It should be 
noted that the muon scintillation detectors are divided into 2 groups placed 
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under different absorber thickness. This allows to study EAS muons with 2 
different energy thresholds: > 5 GeV (Hall) and > 2.5 GeV (Tunnel). 
These thresholds were estimated experimentally using scintillation telescope 
with a small solid angle. The arrangement of the muon detectors gives the 
possibility to determine the muon lateral distribution up to 60m at > 
5 GeV and up to 90m at E^ > 2.5 GeV. 

Some peculiarities of the GAMMA array geometry make necessary to pay 
attention on the azimuth angle symmetry of data. The surface detectors are 
placed on the hillside and sometimes the difference between Z coordinates 
of the individual detectors is reaching 18m. On the other hand, the position 
of muon detector is very asymmetric in regard to the geometric centre of 
our installation. In this case, the influence of azimuth angle ip could be 
strong. To check this effect we divided the observed showers into four groups 
with <p e [0 - 90]°, [90 - 180]°, [180 - 270]°, [270 - 360]°. The comparison 
of the average muon lateral distributions for each of these groups by ip is 
shown in figured (a, b) for < N e >= 1.28 • 10 6 . It can be visible, that the 
difference of the muon lateral distributions for the different ranges of (p is 
quite negligible. Furthermore, it is necessary to point out that for the small 
distances (less than 15 m for >2.5 GeV and 8 m for >5 GeV) there is a 
significant influence of the punch through of particles (high energy electrons 
and hadrons) and these distances are excluded from the further analysis. 



3 Results 

We have used the experimental data obtained during 3300 hours operation 
time. The number of EAS with N e > 10 5 and with zenith angle 6 < 30° is 
^ 260.000. The effective area for the selection of the EAS axis is ^ 5.000m 2 . 

3.1 Electromagnetic component characteristics 

The main parameters of the EAS electromagnetic component at the given 
observation level are: the total number of the charge particles or shower size 
N e and the lateral shower age parameter S. In order to obtain the EAS 
size spectrum and to investigate its behaviour in the knee region the correct 
estimation of these parameters is needed. 
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Usually in the energy region of 10 5 — 10 7 GeVthe scintillation detectors 
are used to register the EAS charged particles. However it is well known that 
such type of a detector registers not only the EAS charged particles, but also 
electrons generated in the absorber above as well as in the scintillator by 
EAS photons. This contribution to the measured charged particle density 
depends both on the scintillator thickness and the distance of detector from 
the shower axis. 

The comparison of the charge particle densities seen by scintillation de- 
tectors and Geiger counters has been made for different distances from the 



shower core P3[. In the same way, comparisons with other kind of detectors 



have been performed, (spark chamber | 24j , Geiger counter [25f], neon flash 



tube p™, as well as with thin and thick scintillators |0[B8j). But, up to 
now, there are not exact estimation of the photon contribution. In the energy 
range 10 5 — 10 7 GeV all these investigations show both a smooth rise and a 
noticeable decrease of the photon contribution to the measured shower size 
versus the distance to the shower core up to 100 m. 

Recently the methodical experiment on the MAKET-ANI installation on 
the Mt. Aragats has been carried out | 29| . Comparing the densities measured 
by scintillators with different thicknesses (1.0, 1.5 and 5.0cm), up to 100m, 
the ratio K 7 (r) of the measured density p sc {r) to the charged one p ch {r) is 
described by K 7 (r) = p sc {r)/ p ch {r) = (r/R M )~ 0A8 , where r is the distance 
from the shower core and Rm = 120 m is the Moliere radius for the obser- 
vation level of 700 g. cm 2 . According to this result the photon contribution 
to the charged particle density is ^ 17% at 50 m from the shower axis and 
practically disappears at 100 m. 

At the same time, simulations show |30|]|3]J a significant increase of the 
ratio of the photon number to electrons one versus r. Calculations for the 
real conditions of the GAMMA installation have been performed taking into 



account the absorber and scintillator thicknesses |32| . According to this work 
the photon contribution to the measured number of the charged particles is 
practically constant from 5 to 100 m and is K 1 ^ 1.25. The processing of 
the GAMMA data considering these two results, K^{r) = (r / Rm)~°' 1s and 

= 1.25, ||29||[j3~2" | , has been made to study the influence of this factor on 
the measured EAS features. 

In figures 4 the experimental electron lateral distributions are shown for 
three EAS average size intervals < N e >=2.3 *10 5 , 7.3 *10 5 and 22.8 *10 5 as 
estimated for the two cases of the photon contribution. It can be seen that 
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in both cases the electron lateral distribution can be very well approximated 
using the NKG function, but because of the different shape of versus r, 
the average age parameter value < S > is different by about 0.16. 

The comparison between the measured and simulated electron lateral 
distributions is shown in figure 5 for the same values of < N e >. It can be 
seen a good agreement for the correction factor K 1 = 1.25 = const. This 
agreement improves with the increasing of < N e >. 

It has to be noticed that the age parameter only describes the shape of 
the charged particles lateral distribution and its values < S > < 1 do not 
contradict to the fact that the shower is over the longitudinal development 
maximum. Figure 6 shows the dependence of the age parameter < S > 
versus the shower size < N e > from our experimental data processed with 
Ky(r) = {t/Rm)~ 8 and = 1.25 and the results from our simulation and 
the Tien-Shan experimental data p3|. For < N e > > (3 — 5) * 10 5 , exper- 



imental data with K n = 1.25 are in a good agreement with the simulation 
results which is not at all the same in the case of K y (r) = (r/i?^) -0 ' 18 for 
which the disagreement is large. 

For the both variants of processing, the differential size spectra are shown 
in figure 7a for showers in the zenith angle interval 9 ^ 30°, (mean at- 
mospheric depth : 738 g. cm" 2 ). In both cases the knee can be easily ob- 
served. Spectra are approximated in the whole N e interval by formula 



I(N e ) = A * N e 71 (1 + (AT e /A^e)fc)-( 7 i-7 2 )/fc. The spectra differ mainly at 
the beginning and the end of the size interval, (by no more than 20%). Spec- 
tral indexes difference for GAMMA for both cases is A7 = 0.36 ± 0.02. The 
knee regions approximately coincide at N^ nee « (1.8 — 2.0)* 10 6 . The integral 
intensity for the range of N e ^ j\f^nee j g ^ g _|_ g ^ I0 _7 m~ 2 s _1 sr _1 . This 
result is in agreement with the data of other experiments ||34j| ||35|| . 

In order to demonstrate the knee region properties more visually spectra 
are given in figure 7b multiplied by N^ 5 . For the case of K y = 1.25 spectral 
index is less by (0.07 ± 0.02) in the whole N e region. 

In the figure 8, the size spectra for the two zenith angle intervals are 
shown at < 9 >= 14° and < 9 >= 31° to make the qualitative assessment 
of the behaviour of the spectra at different angles. They are practically 
parallel and coincide when data with < 9 >= 31° are shifted. Such behaviour 
corresponds to the shower absorption lenght of A = (230 ± 25) g ■ cmT 2 below 
the knee and (230 ± 40) p ■ cm~ 2 above up to iVe = 5 ■ 10 6 . The permanency 
of the attenuation length is one evidence of the constant charged particle 
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composition of the showers below and above the knee only and says nothing 
about the cause of the knee. 

In figure 9 the spectra are shown with various cuts of the age parameter 
S. Relatively old shower spectrum (S > 0.85) practically has no knee and 
is steeper than all shower spectrum. On the other hand, the young shower 
spectrum (S < 0.75) has the obvious knee and is flatter, than the all shower 
spectrum below the knee and is almost parallel to it above. These spectra 
may be explained with both hypothesises of the knee origin noted in the 
introduction, because of the region of N e > (1 — 2) • 10 6 is transitional. The 
young showers fraction does not arise above the knee because of change of 
the mass composition to heavier one in the model of primary knee. In the 
case of change of the hadron interaction the fraction of the heavy-like events 
will rise among the showers generated by primary protons. 

3.2 Muon component characteristics 

As it was mentioned above, the underground muon detectors are placed in 
the hall and tunnel of the GAMMA underground part. Due to the different 
shieldings of concrete, iron and ground above detectors, muons with two en- 
ergy thresholds ^ 5 GeV (hall) and ^ 2.5 GeV (tunnel) are detected. 
The muon lateral distributions for given thresholds and three < N e > inter- 
vals are shown in figure 10 a. These distributions are limited by the number 
of detectors and their disposition in the hall and tunnel. Moreover, close 
to the shower core the signal in the detectors are strongly affected by the 
contribution of other type of particles due to the punch through effect. So, 
the correct determination of the muon lateral distribution is available at dis- 
tances [8 - 52] m for £ M ^ 5 GeV and [20 - 90] m for ^ 2.5 GeV. In 
the interval of [0 — 50] m, the shapes of lateral distribution are very simi- 
lar for the both muon energy thresholds. This gives the possibility to find 
a transition coefficient from the density of muons with E^ ^ 2.5 GeV to 
the muon density with E^ ^ 5 GeV, which is p(2.5 GeV)/p(5 GeV) ~ 1.35, 
for all < N e > intervals. In figure 10 b, the muon lateral distribution 
with the energy threshold E^ ^ 5 GeV for the hall detectors and recon- 
structed one from the distribution of tunnel detectors are presented. These 
distributions are well approximated by the Hillas function, pSf , p (N e ,r) = 
0.9 (iVyiO 5 ) r -(o-75+o.o6-in(Ar e /io5) exp (_ r / 8 o). This formula describes well the 
experimental data for N e < 2.10 6 . Farther, the dependence of p^ versus 
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N e becomes stronger but the shape of the distribution remains the same. 



The approximations of the Tien-Shan data |37|], for the same energy 
threshold ^ 5 GeV and the same observation level 700 g ■ cm~ 2 are 
presented by p M (iV e ,r) = 0.95 (iV e /10 5 ) - 8 r" a75 exp(-r/80). This formula 
describes the GAMMA experimental data noticeably worse, especially for 
N e > 10 6 . 

As it was mentioned above, the lateral distribution of muons in our 
experiment, (E^ ^ 5 GeV), is studied for distances r ^ 52 m from the shower 
core where about 33% of all muons are contained. Obtaining the total muon 
number, iV M ,by integration the muon lateral distribution function, the possi- 
ble extrapolation errors for r > 52 m will give an error in the estimation. 
On the other hand, closer to the shower axis, r < 8 m, other kinds of par- 
ticles give a considerable contribution to muon lateral distribution. For this 
reason, we use the truncated muon number, A^™ n , which is the number of 
muons in the ring 8m < r < 52m. Experimentally, the truncated muon 
number is determined as Nj[ un (r) = YUpT /w^ un (r)) / K , where: 

K is number of detectors in the given interval of r; 

pl x is the experimental muon density in the i th detector and 

u>*™ n (r) is the truncated muon probability distribution, determined by 
our muon lateral distribution approximation. 

It should be noticed, that according to simulations, Nj[ un (>5 GeV) in our 

experiment depends on the primary particle mass. Figure 11 shows N t ^ un 
dependence on N e . Up to N e « 2 ■ 10 6 it can be fitted with the expression 
jytrun ^ jyO- 79 . For larger N e dependence becomes steeper. It should be noted 
that because of the shower size threshold N e > 10 5 the registration of the 
shower with fixed muon size and 100% efficiency is possible at Njf un > 10 3 . 



The differential Nj[ un spectrum is plotted on figure 12. It can be seen 



that there is no obvious knee, but the data do not contradict to change of 
spectral index by (0.1-0.2) at Nj[ un > 5 • 10 3 . Contrary to the case of fixed 
N e the showers at fixed are enriched by the showers generated by heavy 
primaries. If the proton spectrum has a knee only then the A M spectrum will 
have feebly marked knee. At energies above the possible knee position of 
the primary iron spectrum, spectrum will obtain the final index different 
by A7 w 0.5 from present one. It is necessary to expand our measurements 
to this energy region. 

Situation in the case of the change of hadron interactions for sufficientely 
high energies depends on details of the secondary spectra but the difference 
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of indexes will be not bigger. 

In order to interpret our experimental results the instrumental response 
has to be taken into account. For this purpose the detector simulation 
program ARES (ARagats Events Simulation) based on the GEANT 
package. |39| has been developed for the GAMMA installation. The data 
from CORSIKA simulated EAS are used as input for the ARES code. 

To take into account the detector response to EAS muon component 
all secondary particles of EAS at the ground level are passed through the 
absorber and muon detectors. For each shower the deposited energy and 
the muon number in the individual detector are obtained (the methodi- 
cal procedure of muon number estimation is described in |]40|||41|| ). The 
muon lateral distribution are derived using the CORSIKA simulation data in 
the primary energy range of 3*10 5 -10 GeV for the different primary groups. 
Figure 13 (a, b) shows the muon lateral distributions in the shower size range 
1.78-10 5 <N e <3.16-10 5 in the case of the normal mixed composition for the 
tunnel and hall detectors. The experimental distribution and the COR- 
SIKA/ ARES simulation results (with the corresponding muon thresholds) 
are in a good agreement. 



4 Conclusions 

The study of the EAS electrons and muons (E^ > 5 GeV and > 2.5 GeV) 
by the GAMMA array gives reliable information about EAS characteristics 
in the shower size range of 10 5 -10 7 . There are not noticeable contradictions 
of our data with other experiments. Using the contribution of the EAS 
gamma-quanta equal = 1.25 = const an agreement between CORSIKA 
simulation and experimental data is obtained. 

Permanency of the age parameter S at N e > 2-10 6 and steepening oiNf[ un 
depending on N e at the same range by N e conform to the more rapid shower 
development. The data in the shower size interval N^ nee -r 10 • N^ nee are not 
decisive for the problem of the knee origin and enlargement of the study range 
is necessary for the complex installations to make sure conclusion about the 
primary composition and spectra. With the object to extend a possibilities of 
the GAMMA installation it is expected to decrease the registration threshold 
by N e and to enlarge an effective area of the shower selection with N e > 10 6 . 

The present paper is based on the ANI collaboration data bank and ex- 
press the point of view of the given group of co-authors. 
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5 Figures 

Figure 1. Layout of the GAMMA installation 

Figure 2. Schematic view of the muon underground detector 

Figure 3. Muon lateral distributions at the different azimuth angles tp: 

a) hall, b) tunnel 

Figure 4. Electron lateral distributions for two cases of the photon con- 
tribution 

a) K 7 (r) = (r/R M )-°- 18 , b) K 7 (r) = 1.25 = const 

Figure 5. Electron lateral distribution in comparison with the CORSIKA 
simulation for two cases of the photon contribution 

Figure 6. The average age parameter < S > versus shower size < N e > 
at two cases of the photon contribution in comparison with CORSIKA sim- 
ulation and Tien-Shan data 

Figure 7. Differential size spectra at two cases of the photon contribution 
Figure 8. Differential size spectra at different zenith angles 
Figure 9. Differential size spectra at different interval by age parameter 
Figure 10. Muon lateral distribution at different shower sizes in 
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a) hall, b) tunnel 

Figure 11. Average truncated muon number < iV* un > versus number of 
electrons < N e > 

Figure 12. Differential spectrum of truncated muon number 

Figure 13. Muon lateral distributions in comparison with CORSIKA and 
ARES simulation results: 

a) hall, b) tunnel 
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